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We present measurements of elliptic flow and event-by-event fluctuations estab- 
lished by the PHOBOS experiment. Elliptic flow scaled by participant eccentricity 
is found to be similar for both systems when collisions with the same number 
of participants or the same particle area density are compared. The agreement of 
elliptic flow between Au+Au and Cu+Cu collisions provides evidence that the mat- 
ter is created in the initial stage of relativistic heavy ion collisions with transverse 
granularity similar to that of the participant nucleons. The event- by-event fluctu- 
ation results reveal that the initial collision geometry is translated into the final 
state azimuthal particle distribution, leading to an event-by-event proportionality 
between the observed elliptic flow and initial eccentricity. 
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Figure 1. Panel a) V2 versus N par t for Au+Au and Cu+Cu collisions at y/s NN = 22.4, 
62.4 and 200 GeV. Panel b) the dependence of V2 on beam energy for mid-central Au+Au 
collisions. We add to this plot PHOBOS published data for Au+Au at ^/s NN = 19.6, 
62.4, 130 and 200 GeV. The bars in the plots represent statistical errors. 



1. Introduction 

Elliptic flow has been studied extensively in nucleus- nucleus collisions at 
SPS and RHIC as a function of pseudorapidity, centrality, transverse mo- 
mentum and center-of-mass energy 1,2 . One of the most striking observa- 
tions at RHIC is a strong event anisotropy in non-central collisions 3 , which 
is generated through the elliptically deformed overlap region of the collid- 
ing nuclei, resulting in an eccentric distribution of matter and anisotropic 
pressure gradients in the early stages of the expansion 4 . From a micro- 
scopic point of view this strong collective anisotropy is best described under 
the assumption of extremely strong rescattering 5 , strong enough in fact to 
reach the limit of continuum dynamics. To achieve such a strong conver- 
sion of anisotropies from coordinate to momentum space, rescattering has 
to be strong at very early times and local thermalization has to occur while 
the geometric deformation of the source is still large 6 . In this work, the 
comparison of the data from Cu+Cu and Au+Au collisions measured by 
PHOBOS experiment at RHIC provides new information on the interplay 
between initial collision geometry (initial eccentricity) and initial particle 
density in determining the observed final state flow pattern. Studies from 
PHOBOS have pointed out the importance of fluctuations in the initial- 
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state geometry for understanding the large Cu+Cu v2 value. 

2. Elliptic Flow 

PHOBOS has measured elliptic flow as a function of pseudorapidity, central- 
ity, transverse momentum, center-of-mass energy 7 ' 8 ' 9 and, recently, nuclear 
species 10 . In particular, the measurements of elliptic flow as a function of 
centrality provide information on how the azimuthal anisotropy of the initial 
collision region drives the azimuthal anisotropy in particle production. In 
figure la, we show the centrality dependence of V2 at midrapidity (|ry| < 1) 
for Cu+Cu at ^/s NN = 22.4, 62.4 and 200 GeV collision energies, as obtained 
from our hit-based and track-based analysis methods 8,1 °. A substantial 
flow signal is measured in Cu+Cu at three energies, even for the most 
central events. The strength of Cu+Cu v 2 signal is surprising in light of 
expectations that the smaller system size would result in a much smaller 
flow signal . Figure, lb shows the dependence of v 2 on beam energy 12 . 
We add to this plot PHOBOS published data for Au+Au at ^/s^ = 19.6, 
62.4, 130 and 200 GeV 8 . At low fixed target energies (^/s^~ 3 GeV), 
particle production is enhanced in the direction orthogonal to the reaction 
plane, and v 2 is negative. This is due to the effect that the spectator parts of 
the nuclei block the matter in the direction of the reaction plane. At higher 
center of mass energies, these spectator components move away sufficiently 
quickly, and therefore particle production is enhanced in the reaction plane, 
leading to v 2 > 0. This phenomenon is expected in hydrodynamic scenar- 
ios in which the large pressure gradients within the reaction plane drive a 
stronger expansion. However, the most important observation is that, up 
to the highest center of mass energies at RHIC, the observed asymmetry 
v 2 continues to grow. It would be interesting to see whether this tendency 
is confirmed by the LHC data. 

3. Initial Eccentricity 

In order to distinguish collision dynamics from purely geometrical effects, it 
has been suggested that the measured v 2 should be scaled by the eccentric- 
ity of the nuclear overlap area 13 . The PHOBOS collaboration has shown 
that for small systems or small transverse overlap regions, event-by-event 
fluctuations in the shape of the initial collision region affect the elliptic flow. 
Monte Carlo Glauber (MCG) studies have shown that the fluctuations in 
the nuclcon positions frequently create a situation where the minor axis 
of the overlap ellipse of the participant nucleons is not aligned with the 
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- Figure 2. Panel a) V2/(e P art) versus N part for 
Au+Au and Cu+Cu collisions at Ws NN = 200 
■ - GeV. Panel b) shows V2/(e P art} as function 
m ~_ of mid-rapidity (|r)| < 1) particle area density 

• Au-Au200Gev ~ l/(S)(dN/dy) for Cu+Cu and Au+Au collisions, 
o Cu-Cu 200 Gev - The bars in the plots represent the statistical 
■ au-au 62.4 Gev - errors. Panel c) correlation plot of midrapidity 
□ Cn-Cu 62.4 Gev - particle area density l/(S)(dN/dy) as function of 
300 Npart calculated in Glauber MC for Au+Au and 

Cu+Cu collisions at ,/s^7 = 62.4 and 200 GeV. 



impact parameter vector. To account for this effect, PHOBOS has intro- 
duced the participant eccentricity defined as 10 : £ par t = Ji+J^^^ , 
where a xy = (xy) — (x)(y) is the covariance. This definition accounts 
for the nucleon fluctuations by quantifying the eccentricity event-by-event 
with respect to the overlap region of the participant nucleons. For com- 
parison of Au+Au and Cu+Cu collisions at several energies, figure 2a and 
2b show V2/(e pa rt) as a function of A part and midrapidity particle area 
density, l/(S)(dN/dy), respectively. Figure 2c shows that at given energy, 
62.4 or 200 GeV, the Au+Au and Cu+Cu collisions selected for the same 
value of -/Vp ar t leads to the same value of midrapidity particle area density 
l/(S)(dN/dy) calculated in MCG. We observe in figure. 2a (2b) that the 
v 2 scaled by e part are similar for both Cu+Cu and Au+Au collisions at the 
same value of iV part (l/(S)(dN/dy) ). It should be noted that in figure 2b 
which has been introduced previously in Ref. 14 , in the y-axis the V2(?j) has 
been converted to V2(y) by scaling the data by factor 0.9 and also in the 
x-axis the dN/dy = 1.15 dN/dry at midrapidity region, \rj\ < 1. This simi- 
larity between Cu+Cu and Au+Au collisions is also observed as a function 
of transverse momentum as well as in a wide pseudorapidity range 15 . 
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Figure 3. Panel a) {v2) and panel b) cr V2 versus Np^rt for Au+Au collisions at y/sj^ = 
200 GcV. Boxes and gray bands show 90% C.L. systematic errors and the error bars 
represent 1-<t statistical errors. The results are for < r\ < 1 for the track-based method 
and |»7| < 1 for hit-based and event-by-event methods. Panel c) cf V2 /(v2) versus A/part 
for Au+Au collisions at v / s^7 w " = 200 GeV. Open squares show <T £pBrt /(e pa rt} calculated 
in a Glauber MC. The bands show 90% C.L. systematics errors. 



4. Elliptic Flow Fluctuations 

The apparent relevance of the participant eccentricity model in unify- 
ing the average elliptic flow results for Cu+Cu and Au+Au collisions leads 
naturally to consideration of the dynamical fluctuations of both the par- 
ticipant eccentricity itself as well as in the elliptic flow signal from data. 
Simulations of the expected dynamical fluctuations in participant eccen- 
tricity as a function of N part were performed using the PHOBOS Monte 
Carlo Glauber based participant eccentricity model. Figure. 3a shows the 
mean, (v 2 ), and the standard deviation, ct V2 , of the elliptic flow parameter 
V2 at midrapidity as a function of the number of participating nuclcons, in 
Au+Au collisions at y/s NN — 200 GeV for 6-45% most central events 17 . 
The results for (V2) are in agreement with the previous PHOBOS v 2 mea- 
surements 8 , which were obtained with the event-plane method for charged 
hadrons within \rj\ < 1. The uncertainties in dN / drj and v 2 (7/), as well as dif- 
ferences between HIJING and the data in these quantities, introduce a large 
uncertainty in the overall scale in the event-by-event analysis due to the 
averaging procedure over the wide pseudorapidity range. The event-plane 
method used in the previous PHOBOS measurements has been proposed 
to be sensitive to the second moment, \J (vf), of elliptic flow 16 . The flue- 
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tuations presented in this work would lead to approximately 10% difference 
between the mean, (v 2 ), and the RMS, \J (v|), of elliptic flow at a fixed 
value of iVp art . Most of the scale errors cancel in the ratio, cr V2 /( v 2), which 
defines "relative flow fluctuations", shown in figure. 3b as a function of 
the number of participating nuclcons 17 . We observe large relative fluctu- 
ations of approximately 40%. MC studies show that the contribution of 
non-flow correlations to the observed elliptic flow fluctuations is less than 
2%. Figure 3c shows cr epart /(e part ) at fixed values of N part obtained in a 
MC Glauber simulation. The 90% confidence level systematic errors are es- 
timated by varying Glauber parameters as discussed in Ref. 10 . A striking 
agreement between the relative fluctuations in the Glauber model partici- 
pant eccentricity predictions and the observed elliptic flow fluctuations is 
seen over the full centrality range under study. The observed agreement 
suggests that the fluctuations of elliptic flow primarily reflect fluctuations 
in the initial state geometry and are not affected strongly by the later stages 
of the collision. 

5. Summary 

We have performed a comprehensive examination of the elliptic flow of 
charged hadrons produced in Cu+Cu and Au+Au collisions at ^/s NN = 
19.6, 22.4, 64.4 and 200 GeV. We also presented the measurements of 
event-by-event fluctuations for Au+Au collisions at ^/s NN = 200 GeV. The 
comparison of the data from Cu+Cu and Au+Au collisions provides new 
information illustrating that the participant eccentricity is the relevant ge- 
ometric quantity for generating the azimuthal asymmetry leading to the 
observed flow. The magnitude of event-by-event fluctuations agree with 
predictions for fluctuations of the initial shape of the collision region based 
on the Glauber model. These results provide qualitatively new information 
on the initial conditions of heavy ion collisions and the subsequent collective 
expansion of the system. 
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